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ABSTRACT 

We investigate the effects of the gravitomagnetic corrections to the usual gravi- 
tational lens quantities for a specific lensing mass distribution modelled after spiral 
galaxies. An exponential disk is embedded into two different spherical halo models 
where disk and haloes parameters are fixed according to the observed mass to light 
ratios, galaxy magnitudes and rotation curves. The general expressions for the lensing 
deflection angle are given also taking into account the orientation of the galaxy disk 
plane with respect to the lens plane. It is found that the gravitomagnetic term changes 
the deflection angle by a typical amount of the order of ten microarcseconds. 
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1 INTRODUCTION 

Gravitational lensing phenomena are a powerful tool to in- 
vestigate the curvature of space -time nearby the source of 
gravitational field. The mathematical formulation of lens- 
ing is well developed. The weak field approximation may 
be used whenever the distance of the light ray to the lens is 
much larger than its Schwarzschild radius (e.g., ~ 10 times). 
This is the situation in several gravitational systems consid- 
ered in astrophysical applications. When this is not the case, 
the mathematical formalism becomes much more compli- 
cated involving transcendental lens equations which cannot 
be easily solved (Bozza et al. 2001). Furthermore, interest- 
ing contributions to the lensing quantities, in the weak field 
approximation, could come out considering higher order cor- 
rections to the lens potential as the gravitomagnetic ones. 
Generally, such contributions are discarded since are consid- 
ered too small (Schneider et al. 1992). In Capozziello et al. 
(1999), these higher order corrections have been estimated 
for point - like deflectors showing that they could give rise 
to non negligible effects. A further generalization was ob- 
tained by one of us (Sereno 2002). By using Fermat's prin- 
ciple and the standard assumptions of gravitational lensing, 
the gravitomagnetic corrections to the time delay function 
and the deflection angle for a geometrically thin lens were 
derived. The effects of the lens angular momentum on the 
propagation of light rays has been considered in literature 
using different approaches. For instance, Dymnikova (1986) 
determined the time delay of a signal due to the graviational 



fleld of a rotating body integrating the null geodesies of the 
Kerr metric. The effects of the rotation of the deflector have 
been evaluated to the lowest order by Glicenstein (1999) 
who used an argument based on Fermat's principle applied 
to the Lense - Thirring metric. These results are limited to 
pointlike deflectors and are difflcult to generalize to more 
complex geometries. Asada & Kasai (2000) used a multipole 
expansion of the deflecting potential to investigate the case 
of extended lenses. On the other hand, Sereno & Cardone 
(2002) applied the formalism* developed in Sereno (2002) 
to the case of spherically symmetric mass distribution with 
angular momentum in the case of a rigid body rotation. 

The aim of this paper is to take into account the con- 
tribution of gravitomagnetic corrections in the case of spiral 
galaxies modelled according to two realistic mass profiles. 
We study how this contribution affects the lensing deflec- 
tion angle and estimate its detectable range. This approach 
reveals particularly interesting, especially if related to the 
capabilities of the present 10 m class telescopes (like VLT, 
Keck, Subaru), or those ones which are going to be avail- 
able in the next future (like the 100 m OWL telescope) and 
the space based interferometric mission, such as SIM and 
GAIA. The last decade has seen a notably effort in applying 



* This approach has been now further generalized to metric theo- 
ries of gravity (Sereno 2003) to include also post - post - newtonian 
corrections in the computation of the deflection angle. However, 
in this paper, we will not consider these deviations from the stan- 
dard general relativity theory. 
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gravitational lensing as a tool to investigate the dynamical 
structure and morphology of galaxies. In particular, a main 
problem in astrophysics concerns the nature and the distri- 
bution of dark matter in the galactic components, that is the 
halo, the disk or the bulge. Paczyiiski (1986) realized that 
this problem could be faced by gravitational microlensing. 
His predictions have been actually confirmed by the detec- 
tion of several microlensing events towards the Magellanic 
Clouds (Lasserre et al. 2000; Alcock et al. 2001), the bulge of 
our Galaxy (Udalski ct al. 1994) and the Andromeda Galaxy 
(Auriere et al. 2001; Calchi Novati et al. 2002). These can 
be considered as an evidence of the presence of dark (likely 
baryonic) components in the form of massive astrophysical 
compact halo objects (MACHOs). Anyway, MACHOs are 
not the only candidates for baryonic dark matter in the halo. 
It has been proposed, in fact, that an appreciable fraction 
could be in the form of diffuse objects as self gravitating 
gas clouds, with masses of the order of 10~^ A4q and radii 
7? 10 AU (Pfenniger et al. 1994; Kcrins et al. 2002). The 
presence of a dark halo is implied by the flat rotation curves 
in spiral galaxies (Sofue & Rubin 2001). It has been sug- 
gested that the rotation curve of such galaxies follows the 
so - called universal rotation curve that can be characterized 
by a single free parameter, namely the total luminosity (or 
mass) of the disk (Persic et al. 1996a). The profile of these 
curves is given by the sum of an exponential disk and a 
spherical halo with a flat density core. However, this does 
not tell us what is the dark halo mass density, but it allows 
to recover some useful scaling relations to express the halo 
parameters as functions of those of the disk. 

This paper is organized as follows. In Sect. 2, we de- 
scribe the models of spiral galaxies which wc consider and 
the scaling laws that we will use to determine their pa- 
rameters. In Sect. 3, we derive general expressions for the 
deflection angle of a spiral galaxy taking into account the 
gravitomagrictic corrections and the orientation of the disk 
plane with respect to the lens plane. Results arc discussed 
in Sect. 4, while conclusions are presented in Sect. 5. 



2 MODELLING SPIRAL GALAXIES 

Following the standard approach, wc will consider spiral 
galaxies as built by two main components : a disk and a dark 
halo. The disk, essentially, accounts for the optical proper- 
ties of the galaxy. Disk parameters may be estimated by 
fitting the data taken in different photometric bands for a 
given model of the surface brightness distribution. It is usual 
to describe the disk surface density with the exponential 
profile^ (Freeman 1970) : 



E(i?) = Soe 



-R/Rd 



(1) 



Two parameters arc needed to characterize the model : the 
central surface density Eo and the scale length radius Rd. 
The latter is given by fitting the optical data, while So is 



^ We assume the disk is idealised as being infinitely thin. This is 
not a serious shortcoming of our analysis since the finite thickness 
of the disk does not change the main results. 



related to the observable central surface brightness no = 
— 2.51og/o as : 



So = T 7o , 



(2) 



being T the disk mass - to - light ratio (which can also depend 
on the radius). Equation (1) does not tell anything about the 
vertical density distribution, but this further quantity is not 
needed for our aims. The total disk mass is Md ~ 2-k Eo R^. 
Another important quantity is the rotation curve, i.e. the 
circular velocity of stars in the disk potential as function of 
the radius R in the disk plane. For an exponential thin disk, 
it is (Freeman 1970; Binney & Tremaine 1987) : 



Vc,disk{R) = 



2GMd ^ 

Rd 



[Io{y)Ko{y) - h{y)K^{y)] 



(3) 



with y = R/2Rd; In(y),Kn(y) are Bcssel's functions of, re- 
spectively, first and second type of order n. The rotation 
curve increases until it reaches a maximum value and then 
decreases with a keplerian trend. 

It is well known, however, that observed rotation curves 
of spiral galaxies remains fiat well beyond the end of visible 
disk (see, e.g., Sofue & Rubin 2001). Even if other theoreti- 
cal solutions have been suggested (see, e.g., Milgrom 1983), 
the flatness of rotation curves may be easily explained by in- 
troducing a second non - visible component in spiral galaxies 
models : the dark halo. Deriving the shape and mass distri- 
bution of such a component from the analysis of the rotation 
curves is a highly degenerate problem since there are sev- 
eral halo models which flt equally well the observed rotation 
curve of a given spiral galaxy. Some hints on the halo struc- 
ture come from N - body simulations of galaxy formation in 
hierarchical ACDM scenarios. However the results are still 
not unique due to problems of resolution and differences in 
the input basic properties of dark matter particles. 

Among the many proposed distributions, we will con- 
sider here the so - called Burkert - Borriello - Salucci model 
(hereafter BBS; Burkert 1995; Borriello & Salucci 2001; 
Salucci & Borriello 2002), proposed to empirically account 
for the rotation curves of four dark matter dominated galax- 
ies and successfully applied to a sample of 17 dwarf and LSB 
galaxies (Kravtsov et al. 1998). Its density profile is : 



PBBs{r) = 



(4) 



(r+ n,)(r2+r,^) ' 

where po and ro are free parameters which represent the 
central dark matter density and the scale radius. 

Since we adopt spherical symmetry, the halo contribu- 
tion to the circular velocity is : 



Vc,BBs{r) 



GMBBsir) 



(5) 



being MBBs{r) the mass of the BBS halo enclosed within 
r, which is given by : 



Me 



5(r) = 6.4 PO ro^ {in (l + ^) - tan"^ (^) + 



2 



1 + 



Vro/ 



(6) 



The mass distribution diverges so that we need to introduce 
a cut - off radius in order to get a finite total mass. To this 
aim, we will truncate the halo at the virial radius Vvir defined 
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as the radius within which the mean density is Syir Pcrit, 
with 5vir depending on the chosen cosmological model and 
Pcrit = 3Ho/8nG the critical density of the universe. The 
total mass M^ir is thus obtained by putting r = Vvir in 
Eq.(6). To estimate r-yir, wc have to fix the value of Svir 
choosing a background cosmological model. We adopt the 
currently popular ACDM model with {ilrn,^A) = (0.3,0.7) 
when = 337. 

In principle, there is no obvious reason why disk and 
halo parameters should be correlated. Actually, this turns 
out to be true : a clear correlation exists. By considering 
a sample of 1100 synthetic rotation curves built from 15000 
measurements, Persic et al. (1996a,b) have demonstrated the 
existence of the so - called universal rotation curve (hereafter 
URC), a function which describes (within the measurement 
errors) the rotation curve of all spiral galaxies, modulo a 
scaling factor related to the disk total luminosity. Prom the 
very existence of the URC, it is possible to infer some useful 
scaling relations among halo and disk parameters. Briefly, 
one decomposes the URC as : 



2 _ 2 I 2 

VURC — 'fc,disk + ^c^halo 



(7) 



with Vc,disk as in Eq.(3) and Vc,haio determined by the dark 
halo model. The universality of the URC does imply that 
disk and halo parameters should be related by some scaling 
relations in order to reproduce the URC itself. Modelling the 
halo by the BBS profile, Salucci & Burkcrt (2000) have found 
the following scaling relations among the halo parameters 
and the total disk mass (see also Salucci & Borriello 2002) : 

log po = -23.0 - 0.077 log Ma - 9.98x lO"'' Md *^ , (8) 



log n, = 9.10 + 0.28 log po - 3.49x10^° Po■''^ 



(9) 



with Po in g cm , Ma in units of the solar mass M© and ro 
in kpc. Adding to Eqs.(8, 9) a third relation (Persic et al. 
1996a; Salucci & Borriello 2002) : 



log Jid = 4.96 - 1.17 log Md + 0.070 (log Md)^ 



(10) 



it is possible to estimate all the galaxy parameters from the 
knowledge of the disk total mass. In principle, it should be 
possible to use the disk scale length Rd as an order pa- 
rameter since it is directly measurable. However, it is worth 
to remember that it may vary with the photometric band 
chosen to fit the disk model. On the other hand, the disk 
total mass Md may be indirectly obtained by measuring the 
central surface brightness /io and using Eq.(2) since T is in- 
dependently estimated. It is worthwhile to stress that there 
are different possibilities to get an estimate of T. We only 
remember here the correlations between T and the color 
B-R (Bell & de Jong 2001) or between T and Rd (Graham 
2002). 

The BBS model is not the only density profile able to 
reproduce the observed URC. Actually, in their original pa- 
per, Persic et al. (1996a) proposed a different approach to 
the URC by modeling directly the halo contribution to the 
rotation curve. They used : 



(11) 



with Vopt the circular velocity at the optical radius Ropt = 
3.2Rd, X = r/Ropt, a the halo core radius in units of Ropt 
and: 



. 1.1 G Md 

= 2 5 

Kpt Rd 



(12) 



If wc assume spherical symmetry, the mass distribution cor- 
responding to the halo rotation curve in Eq.(ll) is : 



Mpss(x) = 



Ropt'^opt 

G 



(l-/3h)(l + a^) 



x^ + 



(13) 



which is once again divergent. As before, we will truncate 
the halo at the virial radius ryir- The mass inside Vvir is thus 
the total mass of the halo. We will refer, in the following, to 
this model as the PSS model. The density profile turns out 
to be : 



■opt 



(14) 



The URC is then obtained by using the following set of 
scaling relations (Persic et al., 1996a,b) : 



Md 
2nRl 



7.8x10* (l + O.G log-^) 



= 7.8x10* [1 - 0.4 {Mi + 21.9)f , (15) 

I3h = 0.72 + 0.44 log ^ = 0.72 - 0.176 {Mi + 21.9) , (16) 
if* 



a = 1.5 ('^)° ' = 1.5xl0-°-°^« (M,+2i.9) 



(17) 



where Li is the total luminosity (in solar units) in the / 
band, Mr the absolute total magnitude and — 2.51ogL* = 
—21.9. We will use Mi as an order parameter since it can 
be measured from observations; for each value of Mi, we 
first use Eqs.(lO) and (15) to get the disk parameters. Then 
Eqs.(12), (16) and (17) will allow us to estimate the halo 
parameters. So the spiral galaxy model is fully characterized. 

The two halo models described here will be used in the 
following section to estimate whether the gravitomagnetic 
corrections may be directly observed in realistic cases. The 
scaling relations quoted above will help us to estimate the 
effects as a function of only one ordering parameter: the disk 
total mass Md in the BBS case and the disk absolute total 
magnitude Mi in the PSS model. 



3 LENSING DEFLECTION ANGLE 

Following Sereno (2002), the deflection angle to the order 



c reads : 



■I' 



1-2 



vi{a 



(18) 



being vi{^') the weighted average, along the line of sight, of 
the component of the velocity orthogonal to the lens plane, 
i.e. : 



and: 



/ {v{^l)■e^n) P{il)dl 



(19) 



(20) 
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In Eqs. (18), (19) and (20), we are using a coordinate system 
centred on the galaxy lens with the axes (^1,^2) in the lens 
plane which is ortogonal to the direction of the incoming 
light ray (coincident with the line of sight), while I is the 
component along the line of sight. Eqs. (18) and (19) show 
that, in the thin lens approximation, motions in the plane 
of the lens does not give any contribution to the deflection 
angle. 

Due to tidal interactions with other galaxies during its 
formation process, a spiral galaxy has a non-zero angular 
momentum which gives rise to corrections to the lensing de- 
flection angle. The effects of the lens angular momentum on 
the propagation of light rays has been considered by Sereno 
& Cardone (2002). They applied the formalism developed in 
Sereno (2002) to the case of a rotating spherically symmet- 
ric lens in the approximation of rigid rotation, i.e. assuming 
that the angular velocity does not depend on coordinates. 
While useful to estimate the gravitomagnetic corrections for 
lensing by stars, the results in Sereno & Cardone (2002) have 
to be generalized to the case of non costant angular velocity 
to be applied to spiral galaxies. 

To this aim, let us consider a spherically symmetric lens 
that rotates about an arbitrary axis, fj, passing through its 
centre (i.e. a main axis of inertia). To specify the orienta- 
tion of the rotation axis, we need two Euler's angles : /3 is 
the angle between the line of sight / and the line of nodes 
deflned at the intersection of the I ^1 plane and the equato- 
rial plane (i.e., the plane orthogonal to the rotation axis and 
containing the lens centre); 7 is the angle between fj and the 
^2 - axis. Using the axial symmetry about the rotation axis, 
one easily gets (Sereno & Cardone 2002) : 



V • ei„(^i, ^2, = -ijj(R) [Ci COS7 + ^2 cos /3 sin 7] 



(21) 



where w(i?) is the modulus of the angular velocity at a dis- 
tance R = {Ri + from the rotation axis, while .Ri 
(that, given the spherical symmetry of the system, can be 
taken along the line of nodes) and R2 are the axes on the 
equatorial plane. It is : 



Ri = I cos /3 -h sin /3 , 

R2 = — Z sin /3 cos 7-1-^1 cos (3 cos 7 -|- ^2 sin 7 . 



(22) 

(23) 



For a self - gravitating spherically symmetric system, both 
the energy and the angular momentum are integrals of mo- 
tion and the ordered motion of the stars is on the plane 
orthogonal to the rotation axis. Thus the angular velocity 
reads : 



uj{R) = 



Vc{R) _ y/K,diskiR)+-"c,haloiR) 



(24) 



7? R 
Note that the circular velocity entering Eq.(24) is the total 
one, i.e. the sum of the contributions due to both the lumi- 
nous disk and the dark halo. Inserting Eqs. (24) into Eq.(21) 
and the result into Eq.(18), it is straightforward to get the 
following general expressions for the two components of the 
deflection angle : 



ai{^,0) = af^ + ^ (7iaCOS7-|-/i6COS/9sin7) 



8G 



02 i^,6) = af^ + ^ {ha cos 7 -I- hb cos /3 sin 7) 



8G 



(25) 
(26) 



where we have introduced polar coordinates (^, 6) in the lens 
plane. Here, af^ are the two components of the classical 



deflection angle which, for axisymmetric systems, turn out 
to be : 



.(0) 



-inG 



cose I E(e')?'dC' 



^ sine I m')^'d^' 



In Eq.(25) it is: 



(27) 
(28) 



f 

JO 



(C cos 61 -5' cose') cos e' 



e-f e2-2^ecos(0-0') 



- {u;){^',e')de' ,(29) 



poo 

^0 



L 



I (^ cos g -C' cos sing' , , , 

while, in Eq.(26), ha (hb) is defined as ha (hb) by substi- 
tuting the cos 6' with sin 9' in the brackets on the right hand 
side. The average angular velocity^ has been defined as : 



vc{i,e,i)p{^,i) 



dl 



(31) 



having used the axial symmetry of the halo models described 
in the previous section. Note that Vc depends on (^, 6, 1) 
through R because of Eqs. (22) and (23). 

The gravitomanetic corrections we are investigating are 
the terms Iia,b and /2a, 6 that are of order Vc/c. It is now 
straightforward to compute the defiection angle for a spiral 
galaxy considered as a two component system, i.e. an expo- 
nential disk and a dark halo (modelled with the BBS or the 
PSS profile). The results for the two different modelling are 
presented and discussed in the next section. 



4 RESULTS 

To study the detectability of the gravitomagnetic correc- 
tions, we have to determine the images positions in order to 
investigate whether they are significantly modified by these 
higher order terms. This could be a difficult task since the 
lens equations are highly nonlinear and the results depend 
also on the source position. However, the images positions 
are related to the deflection angle so that, if this latter is 
not modified significantly by gravitomagnetic corrections, 
we have no detectable effects on the images angular separa- 
tion which is the easiest quantity to measure. We can thus 
investigate whether the higher order terms introduce appre- 
ciable corrections to the total deflection angle. To this aim, 
we introduce the quantities : 

A,{^,e;Ma) = WO x MU) - ^ ^ = 1,3,(32) 

i Note that, if wc assume solid body rotation, u) is constant and 
may be taken out from the integrals in Eqs. (29), (30). In this case, 
Eqs. (25) and (26) reduce to Eqs. (8) and (9) of Sereno & Cardone 
(2002). 
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Figure 1. log|A| vs logM^; for the BBS model. The curves are 
referred to ^/r„i^ = 0.5, 1.0, 1.5 from top to bottom, while is 
fixed to 7r/4. 

A(C,^;M.)^100x ^^(i.^)-^ (33) 

being a — (af + ai)^^^ the modulus of the deflection an- 
gle. Note that this quantity is now function of 6') and 
not only of ^ since the gravitomagnetic corrections break 
down the radial symmetry of the deflection angle (Sereno 
2002; Sereno & Cardone 2002). In Eqs.(32) and (33), we 
consider only the dependence on the disk mass Md since the 
scaling relations described in Sect. 2 allow us to fully charac- 
terize the model (i.e. disk and halo parameters) as function 
of this latter quantity. For the PSS model, the ordering pa- 
rameter is the disk magnitude, but there is a one-to-one 
correspondence between Mi and A'ld so that we will use Md 
to compare the results for the BBS and the PSS models. To 
estimate A; and A we have to fix the values of {^,9)- Since 
the deflection angle is roughly proportional to the projected 
mass inside we expect the corrections to be most impor- 
tant for larger values of ^. On the other hand, the deflection 
angle also depends on the angular coordinate 9. We have 
decided to compute At and A for three values of ^, namely 
^/rvir = 0.5, 1.0, 1.5, while we fix S = 7r/4 to speed up the 
calculations^ . Finally, we have also to choose the orientation 
of the disk plane with respect to the lens one. We arbitrarily 
fix /3 = and 7 = 7r/4. It is reasonable to expect that choos- 
ing different values of (/3, 7) does not change significantly our 
results. When computing A^ and A, we have included the 
disk contribution to the classical defiection angle, but we 
have neglected the gravitomagnetic corrections for the disk. 
In fact, we have checked that the classical defiection angle 
of the disk is of the same order of magnitude as the gravito- 
magnetic corrections for the halo at the positions where we 
calculate Ai and A. Hence, the gravitomagnetic corrections 
for the disk will be much smaller than those for the halo and 
can be safely neglected^. We stress, however, that the disk 

§ Note that, for this value of 6, we get : Ai = A2 and A = \/2Ai- 
We do not explicitely calculated Ai and A for other values of 
since it is reasonable to expect that A(5, 8 it/4) is of the same 
order of magnitude as A(^,6 = 7r/4) that is enough for our aims. 
^ Actually, this could not be true for high inclination angles : 
disks almost perpendicular to the plane of the sky may have such 



Figure 2. Same as Fig. 1 but for the PSS model. 



plays a key role in our analysis. First, it is the disk mass 
Md which, through the scaling relations discussed in Sect. 2, 
determines the halo parameters and hence the gravitomag- 
netic corrections. Second, the disk potential contributes to 
the total circular velocity entering Eq.(31). It is thus impor- 
tant to include the disk in our analysis, although we neglect 
its contribution to the gravitomagnetic terms^. In Figs. 1 
and 2, we show A as function of the disk mass Md; a log- 
log plot is used to better illustrate the results**. Despite the 
different halos and scaling laws used, the results for the two 
models are qualitatively similar. This is not an unexpected 
result. Eqs.(25) - (31) show that the gravitomagnetic cor- 
rections depend on the circular velocity Vc{R) of the model 
averaged along the line of sight. On the other hand, the pa- 
rameters of both the BBS and the PSS models have been 
fixed so that the two models reproduce the same URC. As a 
consequence, for a fixed value of Md, the functional form of 
Vc{R) and thus of 9) is the same and hence A has the 

same trend for the two different halo models. The magni- 
tude of the percentual correction is however quite different, 
the PSS one being almost two orders of magnitude higher. 
This result reflects the difference between the BBS and the 
PSS model. Qualitatively, A cx 1/a ~ l/a'°> oc M~^(^) and 
hence the higher is the mass enclosed within ^, the lower is 
the percentage correction. Since, for a fixed Md, the value 
of M{£) for the BBS model is higher than that for the PSS 



a high projected surface mass density that they play a dominant 
role in the appearance of lensed image configurations (Mailer et 
al. 1997; Bartelmann &; Loeb 1998). For these configurations, we 
should take into account also the disk contribution to the grav- 
itomagnetic corrections. However, it is quite unlikely that this 
alters significantly our main results since defiection angles will be 
estimated at J > 0.5 r„ir that is far from the disk edge since 
Tvir >> Rd- For these values of the impact parameter, the disk 
rotational velocity has almost vanished and hence the gravito- 
magnetic corrections due to the disk could be neglected. 
II This also saves computational time without introducing any 
significative systematic error. 

** Note that A(5,9) is always negative for the cases reported 
in Figs.l and 2 so that the plots report log|A|. However, this is 
not a general result depending on the particular choice of the 
orietantion angles (/3,7). 
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halo due to the fact that the BBS mass distribution is more 
concentrated, it turns out that A is larger for the PSS model. 

As a final remark, wc note that the formulae in Sect. 2 
may be used, in principle, whatever is the astrophysical sys- 
tem considered provided that its surface mass distribution 
S(5) is given and the system is spherically symmetric. A typ- 
ical example is a cluster of galaxies which could be modelled 
using the NFW density profile (Navarro, Prenck and White 
1997) and opportune scaling relations. Since the gravitomag- 
rictic corrections roughly scales as Vc/c, one could suppose 
that A should be higher for a cluster of gala:x;ies rather than 
for a single spiral galaxy. A simple scaling argument may 
qualitatively show that this is not true at all. The ratio of 
th gravitomagnctic corrections is roughly : 



J galaxy 



{V/C)gala 



10" 



while for the total deflection angle we get approximately : 
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Thus we obtain : 

^galaxy _ [('^ ) ^ . 

^cluster [("-a<"')Hc;nate 

SO that the gravitomagnetic corrections for a cluster of galax- 
ies are two orders of magnitude lower than for a spiral 
galaxy. Some attempts, using Eqs.(18)-(31) and the NFW 
density profile, have strenghtened this qualitative result. 



5 CONCLUSIONS 

The detectability of the gravitomagnctic corrections wo have 
evaluated is a hard task. Figs. 1-2 show that, for the BBS 
model, |A| may be as high as 0.002%, while, for the PSS 
model, the corrections ranges up to 0.16%. These values re- 
fer to percentage corrections. It is thus useful to give some 
absolute estimates of the gravitomagnetic corrections. For 
the BBS model we got : 

|a-Q°| ~21^857 Atas for CA-„ir = 0.5 , 

|a — a°| ~ 5 ^ 71 fias for ^/r^ir ~ 1.0 , 

\a — a°| ~ 2 7 jias for S^jr^ir = 1.5 , 

while for the PSS model it is : 

|a — q;°| ~ 29 350 p.aa for (^jrmr = 0.5 , 

|a-Q;°| ~ 8^94 /xas for 5/r„i,. = 1.0 , 

|a — a°| ~ 3 43 [las for ^/r„ir = 1.5 , 

having evaluated all these quantities for Q = tt/4 and 
(/3, 7) — (0,7r/4). The images positions depend implicitely 
on the deflection angles so that it is difflcult to estimate the 
variation of the images angular separations A9 due to the 
higher order terms. If one assumes that the correction to 
q(^, 6) leads to a similar variation of A^, then the prospect 
to detect the effects of the gravitomagnctic terms could be 
considered quite good. Actually, the astrometric precision of 
the NASA Space Interferometric Mission (to be launched in 
2009) is estimated to be ~ 4 fms, while the european GAIA 



satellite (scheduled for 2010) will be able to measure stars 
position with an accuracy of ~ 1 j-ias. Both these satellites 
will be thus able to detect tiny deviations in the images po- 
sitions from the values predicted by the standard lensing 
theory. These instrumentations can thus lead these higher 
order terms in the realm of detectability. 

In addition to the rotational velocity, galaxies typically 
have also a peculiar line - of - sight velocity whose effect on 
the deflection angle is easy to estimate. We simply write the 
total deflection angle as the sum of the standard term, a'-'^-' 
and the gravitomagnetic correction, a^^'^ . This latter may 



then be splitted as a 



-ha^ec , being a^fj 



the contribute of the lens angular momentum and cip^c the 
term due to the peculiar motions. Sereno (2002) has shown 
that: 

a°^^ = -2a(°'^ . 

For a receding lens, Vpec < so that the effect of the lens pe- 
culiar motion increases the gravitomagnetic correction thus 
strenghtening the possibility of detection. 

On the other hand, the gravitomagnetic corrections due 
to the dark halo angular momentum could also be lower 
than we have evaluated. A key ingredient in the estimate of 
the higher order terms through Eqs.(29), (30) is the average 
angular velocity (w) that we have computed using Eqs. (24) 
and (31). This approach implicitely assumes that the halo 
objects have only an ordered motion. While this is a correct 
assumption for the disk particles, the high value of the halo 
velocity dispersion suggests that random motions could play 
a significative role in determining the angular momentum 
J. As a test, one could evaluate the spin parameter A = 
J\E\^I'^ IGMI{1 , with E the total energy (Mo et al. 1998), 
and compare it to the results of the simulations of galaxy 
formation (see, e.g., Vitvitska et al. 2002). We find that A is 
overestimated by a factor of order 7 thus suggesting that 
.7 is overestimated too. Should this come to be true, the 
values of the gravitomagnetic corrections we have reported 
should be lowered by a factor of order ~ 1/7. However, it 
is also possible that the angular momentum of the nowaday 
halo is higher than that of the protogalactic halo because of 
the exchange of angular momentum between the halo and 
the collapsed baryons. 

Even if we neglect this question, however, it is worth 
noting that the variations of the deflection angle due to the 
gravitomagnetic corrections are smaller or of the same or- 
der as the ones due to more classical effects. Uncertainties in 
the value of the model parameters may induce errors on the 
value of a'"^ and these latter could be so high that the grav- 
itomagnctic corrections are completely masked. Even if the 
model parameters were known with high accuracy through 
independent techniques, an external perturbing shear (due 
to, e.g., the cluster of galaxies which the main lens belongs 
to or to large scale structure) could mimic the same effect as 
the higher order terms we have considered in this paper. On 
the other hand, this problem may also be reversed. There 
are many multiply imaged QSOs systems where there are 
no evidences of a cluster of galaxies, but nonetheless a small 
external shear is needed to fit the observed images configu- 
ration. Deviations from axisymmetry in the lens galaxy are 
usually claimed to be the origin of this perturbing shear, 
but there are no definitive evidences favouring this hypoth- 
esis. It should be interesting to reconsider these systems to 
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see whether the gravitomagnetic corrections (whose eflects 
could be of the same order as that of a small shear) could 
reconcile model and data without invoking this external per- 
turbation. 

Our results are encouraging since they suggest that the 
gravitomagetic corrections may indeed be detectable. How- 
ever, a further analysis is needed in order to investigate 
how the images position is modified and whether corrections 
emerge also for other lensing observables (such as the ampli- 
fication). Actually, it is especially interesting to consider the 
effect of higher order terms on time delay At between the 
images since this is one of the most promising direct tech- 
nique for the measurement of the Hubble constant (Rcfs- 
dal 1964; Cardone et al. 2002). The time delay At roughly 
scales with the angular separation between the images 
(see, e.g., Schneider et al. 1992). Since the gravitomagnetic 
corrections alter the lensing deflection angle, these higher 
order terms could change and thus At leading to a dif- 
ferent estimation of -ffo- Furthermore, the reconstruction of 
the lensing potential from images configuration can be af- 
fected by the gravitomagnetic corrections since the deflect- 
ing potential tp{^), which is used to evaluate the quantity 
to be fitted, is modified (Sereno 2002). Neglecting gravito- 
magnetic corrections thus introduces a systematic error in 
the determination of the Hubble constant. Investigating in 
detail this issue will be the topic of a forthcoming paper. 
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